Closing the Window on Strongly Interacting Dark Matter 
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Constraints are placed on the spin-independent interaction cross section of dark matter with 
regular matter by refining two methods. First, dark matter-cosmic ray interactions are considered, 
wherein cosmic ray protons collide with dark matter to contribute to the gamma ray sky. This 
constraint is developed using the NFW and Moore dark matter density profiles and new data from 
the Fermi gamma ray space telescope. Second, the Earth capture scenario is considered, wherein 
particles that are captured self-annihilate at Earth's center, thus adding to its internal heat flow. 
The constraint presented here is developed based on analysis of the drift time of dark matter 
particles through Earth, modeled as a core composed of iron and a mantle composed of oxygen with 
linear density gradients between layers. An analysis of the cosmic ray constraint (which rules out 
dark matter-regular matter interaction cross sections greater than its value) shows that it overlaps 
significantly with the Earth drift time constraint (which rules out cross sections smaller than its 
value), closing the window on strongly interacting dark matter particles up to a mass of about 10 17 
GeV when combined with other exclusions. 

PACS numbers: 95.35.+d, 95.30.Cq, 98.70.Sa 
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I. INTRODUCTION 

While the evidence corroborating dark matter's ex- 
istence continues to accumulate [l| @ @, the particle 
nature of dark matter is still unknown. It is impor- 
tant to formulate models of specific candidates, but it 
also is useful to keep the nature generic. As summa- 
rized nicely by Taoso, Bertone, and Masiero 3, dark 
matter at its simplest must be a neutral, cold, weakly 
interacting massive particle that existed in the early uni- 
verse. While a weak interaction strength (corresponding 
to small dark matter-regular matter interaction cross sec- 
tions) is strongly suggested, larger cross sections have not 
been completely ruled out. The parameter space of the 
spin-independent dark matter-nucleon interaction cross 
section a X N and the dark matter mass m x continues to 
be whittled down, however. 

Dark matter's interactions with regular matter can be 
grouped into three categories relative to Earth. In the 
first, dark matter would scatter with particles in the 
Galactic environment, and therefore not interact with 
Earth, resulting in Galactic-scale effects. In the second, 
dark matter would interact weakly enough to make it 
through the Galaxy, through Earth's atmosphere, and 
also through a few kilometers of Earth's surface before 
scattering with regular matter. Underground detectors 
like XenonlOO @, and others like CDMS @ 0, CRESST 
Q , and Edelweiss Q [l(| , have ruled out areas of the dark 
matter mass-cross-section parameter space in this region. 
A third regime lies between the first two, where dark mat- 
ter can pass unaffected through the Galaxy, but interacts 
with Earth's atmosphere or the first bit of Earth's sur- 
face. Indirect detectors like IceCube [111 an d SkyLab [l2| 
[l3| have helped rule out areas of the parameter space in 
this region. 

The findings in Mack, Beacom, and Bertone [l4| 
strongly suggest that dark matter must indeed be weakly 
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FIG. 1. In this log-log plot, the shaded regions of the mass- 
cross-section parameter space shown here are excluded (m x 
is the dark matter mass and g x n is the dark matter-nucleon 
interaction cross section). Other constraints exist, but are 
covered by the shown regions. The values displayed are those 
prior to the analysis of this paper. The astrophysical limit 
from cosmic rays 1 151 ] and the latest underground detector 
limits (XenonlOO are shown in green and orange, respec- 
tively. The middle regions are those from an analysis of data 
from SKYLAB [3 M in red, IceCube O in blue, and the 
Earth heat scenario [1J1 m purple). This paper concerns the 
green (upper-left) and purple (left-center) regions. 



interacting over a wide range of dark matter masses due 
to analysis of this third regime. The upper edge of the 
newly-constrained region from that investigation (purple 
region labeled "Earth Heat" in Fig. [1} lines up with the 
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lower edge of the constraints of the first category (solid 
green line in the same figure) , ruling out a wide swath of 
dark matter interaction strengths over all three regimes. 
However, a more forceful statement towards closing the 
window on strongly interacting dark matter can now be 
made since further research presented here has caused 
these constraint regions to overlap significantly. Note 
that in this paper "strong" interactions refer simply to 
those that are significantly higher in value than "weak" 
interactions. 



II. COSMIC RAY CONSTRAINTS 

The Galactic contribution to the gamma ray sky comes 
primarily from the pion production reaction of cosmic ray 
protons colliding with interstellar hydrogen nuclei (pro- 
tons) . By analogy, collisions between cosmic ray protons 
(p) and dark matter particles (x) would also produce 
gamma rays of similar energy: 



X + P 



X+p + 2-f 



Following Cyburt et al. [15|, the emissivity g 7 (number 
of photons per volume per time) is given by 
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2p(r)^* f 



(1) 



where <£> p is the angle-integrated cosmic ray proton inten- 
sity for all energies above the pion production threshold. 
Note that in this section, p(r) refers to the galactic dis- 
tribution of dark matter. 

This paper uses the Navarro- Frenk- White (NFW) [l6j] 
and Moore [13] [l8| dark matter density profiles for p(r), 
for r in kpc. These profiles are of the form 



p{r) 



Po 



with the following parameters: 

TABLE I. Dark matter density profiles 



Profile 


Q 


P 


7 


r s (kpc) p{R sc ) (GeV/cm 3 ) 


po(GeV/cm 3 ) 


NFW 


1 


3 


1 


20 0.3 


0.259 


Moore 


1.5 


3 


1.5 


28 0.27 


0.0527 



The po term is chosen to ensure that the local density 
(at r — Rsc, the radius of the solar circle) of dark matter 
is 0.3 GeV cm~ 3 for NFW (or 0.27 GeV cm~ 3 for Moore). 



Following Cyburt et al. 



the proton intensity $„ is 



taken to be a constant value in the high dark matter 



mass limit (m x 
For smaller m 



— > oo), so that <£> p 
this scales as: 



11.8 cm" 



*p = $p,oo 



m^o (2m p + m n o ) 
2m x (m p + m^o) 



(2) 



where r ~ 1.7 [151] . Here $ p represents the intensity of 
cosmic ray protons that are energetic enough to produce 
gamma rays. 

The gamma ray intensity due to these dark matter 
cosmic ray interactions, viewed in any given direction, is 
the line integral of q 7 along that line of sight I: 



2n m v p 



p(r)dl. 



(3) 
(4) 



The NFW density profile p(r) is expressed in terms of the 
Galaxy's radial coordinate. Transforming this coordinate 
so that Earth is the origin for a line of sight through the 
Galactic center gives 



Rsc 



21R SC + 1 2 



(5) 



where Rsc = 8.5 kpc, the radius of the Solar Circle (or 
the distance from the Sun to the Galactic Center) (18f . 
The NFW profile in these coordinates is 



P(l') = 



( V72.25-17i'-H' a \ 

V 20 J 
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v / 72.25-17i'+i /2 
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with I' being i/kpe. 

To avoid the divergence of the integral in Eq. 0] at 
the Galactic Center (/ = Rsc), the integral is split into 
three regions, with a central region of radius 0.015 kpc 
treated as having a uniform density p c of dark matter, 
as in Yuksel et al. 18], such that p c — p(8.485 kpc) = 
p(8.515 kpc). The gamma ray intensity is then 
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The upper limit on the integral, 21 kpc (or a distance of 
12.5 kpc from the Galactic Center along the line of sight), 
marks the extent of the the gamma ray region of interest; 
the intensity outside this region is not high enough to be 
significant [19| . 

The diffuse gamma ray emission (DGE) intensity in 
the direction of the Galactic Center is 2.7 x 10 -4 cm~ 2 
sr _1 s _1 when averaged over the four models presented 
by Ackermann et al. [n| based on the observations of 
Fermi-LAT. This value is further refined by accounting 
for contributions to the DGE by cosmic ray interactions 
with Hi, Hii and H2 clouds. On subtracting these con- 
tributions using an average of the models presented by 
Ackermann et al., the maximum intensity attributable to 
dark matter interactions in the direction of the Galactic 
Center is 3xl0 -5 cm~ 2 sr _1 s . However, Ackermann et 
al. also show that pion production is responsible for only 
about two-thirds of the DGE, rather than all of it. This 
gives a final </> 7 value of 2 x 10 -5 cm -2 sr _1 s -1 . Thus, 
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in the high dark matter mass limit, evaluating Eq. [6] for 
the NFW profile gives 



= 4.9 x 10" 29 cm 2 GeV~\ 



(7) 



Any values of a x pj /m x larger than this will result in too 
much gamma ray emission; values lower than this are 
allowed. If we relax the constraint that m x — > oo, then 
from Eq. [2] we have 



1 



m^o (2m p + m,o ) 
2m x (m p + m^o) 



(8) 



where is the cross section in the high dark matter 
mass limit as shown in Eq. [7J This nonlinear behavior 
at low m x values is due to the mass dependence of $ p as 
discussed above. This constraint (including the nonlinear 
behavior) is displayed in Fig. [2] 

The Moore profile, with Earth as the origin for a line 
of sight through the Galactic Center, is 



p{l') 



V72.25-17l'+l' 2 
28 



V72.25-17i'+i' : 
20 



1.5 



with V being Z/kpc. Integrating to a limit of I = 21 kpc 
in three parts (as for the NFW profile) with the same 



value of cf> 7 , we have 



a xN 



9.3 x 1CT 30 cm 2 GeV" 



(9) 



which is a lower limit than that obtained for the NFW 
profile, as expected. 

If instead the full DGE value of 2.7 x 1CT 4 cm" 2 sr" 1 
s _1 were used, the values for a X N/m x in the high dark 
matter mass limit become 6.6 x 10 cm GeV" 1 and 
1.2 x 10" 28 cm 2 GeV- 1 for the NFW and Moore profiles, 
respectively. 



III. EARTH DRIFT TIME CONSTRAINT 

Mack, Beacom, and Bertone [l4[ explored the third 
regime mentioned in Section U by investigating the 
consequences of dark matter's interaction with regular 
matter — namely, dark matter's loss of energy and conse- 
quential gravitational capture by Earth. The dark mat- 
ter would sink to the bottom of the gravitational well, 
Earth's core, where it would annihilate with other dark 
matter particles that experienced the same fate. (In that 
analysis as well as this, dark matter is assumed to be 
its own antiparticle.) If enough dark matter is collected 
through this process the subsequent annihilation would 
produce a heat flow coming from Earth's core that vi- 
olates the measured value. This situation defines the 
lower edge of the exclusion region: any spin-independent 
dark matter-regular matter interaction cross section u x n 
above this line would result in too much heat. (Note that 



in Fig. [T] the bottom line of the constraint from Mack, 



Beacom, and Bertone [14( has been shifted downward 
from the original by a factor of two to reflect a calcula- 
tion^ error in the original work.) 

The upper edge of the exclusion region is where the 
scenario breaks down: the dark matter would interact 
so frequently with the regular matter that it would not 
make it down to the center in a sufficient timescale and 
therefore would not efficiently annihilate. The original 
analysis followed Starkman et al. 20:] , requiring that the 
dark matter be able to drift down to the bottom of the 
gravitational potential well in a time period of 1 Gyr. 
That analysis also assumed an interaction of dark matter 
with iron nuclei. We have refined this calculation and 
show that the new value is very similar to the original. 

In order to more accurately represent the drifting of 
dark matter to Earth's center, we break its interior into 
layers, each with a corresponding size, density, and domi- 
nant target nucleus for interaction. Following Mack, Bea- 
com, and Bertone, to prevent the unlikely consequence 
that any dark matter particle merely coming into con- 
tact with Earth will be captured (and to exclude glanc- 
ing trajectories, for instance), only those particles that 
travel a path length of at least 0.2 R® are considered to 
be captured [14]]. Thus, the effective scattering volume 
for Earth has a radius of 0.99 R®, or roughly R = 6300 
km. 

Four layers of the Earth were considered in order to 
accurately model dark matter interactions in its interior: 
the inner core (i), the outer core (o), the lower mantle 
(1), and the upper mantle (u), as shown in Table 1. This 
analysis treats the core layers as composed entirely of 
iron, and the mantle layers as composed entirely of oxy- 
gen. The density distribution of each layer was treated as 
a linear gradient between the values at the layer bound- 
aries. Density values are approximated from the Prelim- 
inary Reference Earth Model (PREM) (U, and temper- 
ature values are conservative estimates of average layer 
temperatures. Each gradient produced a layer mass con- 
sistent with PREM. 



TABLE II. Earth model setup 



Layer 


Height (km) Density (g/cm s ) 


Temp (K) 


Inner Core (i) 


- 1200 


13.1 - 12.8 


6000 


Outer Core (o) 


1200 - 3500 


12.2 - 9.9 


5000 


Lower Mantle (1) 


3500 - 5700 


5.6 - 4.4 


3500 


Upper Mantle (u) 


5700 - 6300 


4.0 - 3.4 


1500 



For a dark matter particle in equilibrium, with a con- 
stant drift velocity, the gravitational force on it due to 
Earth equals the drag force on it: 



GM(r)m x 



v x AVthnm p v d 



(10) 



where G is the universal gravitational constant, M(r) is 
the mass of the earth interior to the particle at radius 
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r, m x is the mass of the dark matter particle, <j x a is 
the interaction cross section of dark matter with target 
particles of mass number A, vth is the thermal speed of 
the target particles, n is the number density of the target 
particles, rriA is the mass of the target particle, and Vd is 
the drift velocity of the dark matter particle. 

The drift velocity can then be expressed as distance 
traveled over drift time (vd = f/rd), and the mass of the 
target multiplied by the number density can be expressed 
as the density of target particles (nrriA — p{r)), while the 
thermal velocity of the target particles can be expressed 
as that of an ideal gas (i>th = \J 6kT/m,A) with the tem- 
perature being conservatively estimated as an average for 
each layer. Eq. [10] can then be expressed in terms of the 
drift time. 

Note that p in this section refers to the density of the 
target particles in Earth. Since p is a function of radius, 
a small time interval drd in which the particle drifts a 
distance dr is 



dTd=r-^Y-£ 

m x JG 



dr \ M(r) 



dr. 



(11) 



where the density function p(r) for each layer was de- 
rived by constructing a linear gradient between the layer 
boundaries. 

Integrating Eq. [IT] between the layer boundaries us- 
ing the density gradients gives a value for the drift time 
through a specific layer, in terms of a x A/rn x . The in- 
tegral was taken until a radius of 100 m, which marks 
the effective radius of a collected volume of dark mat- 
ter, where the particle can annihilate [l4[ [11] • F° r the 
timescale imposed, the total of all four r values must be 
1 Gyr: 



1 Gyr 



n + T + n 

a x Fe i 



(12) 



(7.66 x 10 12 kg s m~ z ) 



(13) 



+2^(3.76 xlO 12 kgsm- 2 ) 
m x 

This more usefully should be put in terms of a X N/m x , 
which is the interaction cross section of dark matter with 
a single nucleon, rather than a specific target: 



t x n 



A 2 



<j x nA 2 



(t x nA a 



PA 
PN, 

m x rriA 
m x + rriA 

m x + mjv 
m x + rriA 



m x + tojv 
m x mN 



(14) 
(15) 
(16) 



where p is the reduced mass of the dark matter particle 
and the target, which is either a nucleus (A) or a nu- 
cleon (N), and A is the mass number of the target parti- 
cle. The A 2 term indicates that at these low momentum 
transfers, the target nucleus is not resolved, and it is as- 
sumed that the dark matter particle couples coherently 



to the net "charge" or the number of nucleons [14j . Note 
that in the high dark matter mass limit (m x >> tua), 
the enhancement term reduces to A 4 . 

From Eq. 1131 and Eq. 1161 for m x » ttla using A = 56 
for Fe and A = 16 for O this becomes 



a xN 



7.4 x 10" 27 cm 2 GeY -1 . 



(17) 



Values of <j x n /m x lower than this result but greater than 
the bottom edge of the Earth heat constraint in Ref. [l4| 
are excluded; values greater than this result are allowed. 
The behavior of <j x n when the high dark matter mass 
limit is relaxed is shown in Fig. [2] The numerical value 
varies insignificantly from the original result of 7.8 x 10~ 27 
cm 2 m x /GeV, but the behavior at low masses reflects 
the addition of oxygen as a target of interest. This very 
slightly alters the curvature of the limit around m x = 16 
GeV, which reflects the fact that the most efficient scat- 
tering occurs when the dark matter particle mass equals 
the target nucleus mass. The insignificant variance of 
this result from the previous version shows the robust- 
ness of the calculation. Note that this value is more than 
two orders of magnitude above the NFW profile limit ob- 
tained for the cosmic ray-dark matter interactions in the 
previous section. 
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FIG. 2. In this log-log plot, updated exclusion regions as 
determined in this paper are shown in green (upper-left) and 
purple (left-center) with solid border lines. Shaded regions are 
excluded. The cosmic ray constraint shown here (upper left) 
is based on the NFW profile, but it accounts for contributions 
to the DGE by interstellar gas clouds and the pion production 
contribution. Note that using the Moore profile would lower 
this constraint by a factor of ~5, increasing the overlap. The 
cosmic ray constraint not only overlaps with the Earth drift 
time constraint, but also rules out large cross sections for 
dark matter masses up to ~10 17 GeV when combined with 
the SkyLab, IceCube, and underground detector constraints. 
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IV. CONCLUSIONS 

For the Earth capture and drift time analysis, a con- 
straint of cr x jv = 7.4 x 1CP 27 cm 2 m x /GeV was ob- 
tained (in the high dark matter mass limit) as the upper 
limit of all cross sections that would lead to Earth over- 
heating due to dark matter accumulation and annihila- 
tion. By comparison, the original Earth drift time con- 
straint obtained by Mack, Beacom, and Bertone [l4| was 
7.8 x 10 -27 cm 2 m x /GeV in the high dark matter mass 
limit. The insignificant variance in the results shows the 
robustness of the situation and calculation. 

Based on Fermi-LAT observations of the gamma ray 
sky, taking into account the contributions to the DGE 
by cosmic ray interactions with Hi, Hn, and H 2 and the 
fraction of the gamma rays produced by pion decay, a 
constraint of <j x m = 4.9 x 10~ 29 cm 2 m x /GeV (in the 
high dark matter mass limit) was obtained for the NFW 
profile; cross sections larger than this value are excluded. 
The Moore profile provides a stronger overlap, as ex- 
pected, yielding a value of 9.3 x 10~ 30 cm 2 m x /GeV for 
the same high dark matter mass limit. By comparison, 
the original cosmic ray interaction constraint obtained by 
Cyburt et al. [H was 4.6 x 1(T 27 cm 2 m x /GeV in the 
high dark matter mass limit, using an isothermal dark 
matter profile. 

As shown in Fig.[5J the NFW constraint (which is more 
conservative than the Moore constraint) is more than two 



orders of magnitude smaller than the upper edge of the 
Earth drift time constraint. This is a strong statement 
excluding all cross sections in that region, definitively 
ruling out strongly interacting dark matter up to a mass 
of 10 10 GeV (since the Earth capture scenario cannot be 
applied above this mass [14j ) . When combined with the 
SkyLab, IceCube, and underground detector constraints, 
these results serve to rule out strongly interacting dark 
matter up to a mass of ~ 10 17 GeV. Note that if the more 
conservative approach of using the full value of the DGE 
were taken for comparison, the cosmic ray limit still over- 
laps the Earth drift time limit by an order of magnitude, 
ruling out strongly interacting dark matter up to a mass 
of 10 10 GeV. The results in this paper show more defini- 
tively that the dark matter-nucleon spin-independent in- 
teraction cross section must indeed be weak for a very 
large range of dark matter masses. Further analysis can 
be done looking away from the Galactic Center, where 
the gamma ray emission falls much more sharply than 
the dark matter distribution. 
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